Determination of the phase of an electromagnetic field 
via incoherent detection of fluorescence 
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We show that the phase of a field can be determined by incoherent detection of the population 
of one state of a two-level system if the Rabi frequency is comparable to the Bohr frequency so 
that the rotating wave approximation is inappropriate. This implies that a process employing the 
measurement of population is not a square-law detector in this limit. We discuss how the sensitivity 
of the degree of excitation to the phase of the field may pose severe constraints on precise rotations 
of quantum bits involving low-frequency transitions. We present a scheme for observing this effect 
in an atomic beam, despite the spread in the interaction time. 
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It is well known that the amplitude of an atomic 
state is necessarily complex. Whenever a measurement is 
made, the square of the absolute value of the amplitude is 
the quantity we generally measure. The electric or mag- 
netic field generated by an oscillator, on the other hand, 
is real, composed of the sum of two complex components. 
In describing semiclassically the atom-field interaction in- 
volving such a field, one often side-steps this difference 
by making the so-called the rotating wave approximation 
(RWA) , under which only one of the two complex com- 
ponents is kept, and the counter-rotating part is ignored. 
Under this approximation, an atom interacting with a 
field enables one to measure only the intensity, and not 
the phase of the driving field. This is the reason why 
most detectors are so-called square-law detectors. 

In this article, we show how a single atom by itself can 
detect the absolute phase of a driving field, by making use 
of the interference between the co- and counter-rotating 
parts of the excitation, while the Rabi frequency is not 
negligible compared to the transition frequency. This de- 
tection is performed by measuring incoherently the pop- 
ulation of either of the two states of a two level atom. 
This implies that a process employing the measurement 
of population is not a square-law detector in this limit. 
We discuss how the sensitivity of the degree of excita- 
tion to the phase of the field may enable phase telepor- 
tation using a pair of entangled atoms, but poses severe 
constraints on precise rotations of quantum bits involv- 
ing low-frequency transitions. We also present a scheme 
for observing this effect in an atomic beam, despite the 
spread in the interaction time. 

We consider an ideal two-level system where a ground 
state |0) is coupled to a higher energy state |1). We also 
assume that the «-> 1 transitions are magnetic dipolar, 
with a transition frequency ui, and the magnetic field is of 
the form B = B$ cos(tot + </>). We now summarize briefly 
two-level dynamics without the RWA. In the dipole ap- 
proximation, the Hamiltonian can be written as: 



H = e(a -a z )/2 + g(t)a x (1) 

where g(t) = —go [exp(iujt + i<f>) + c.c] /2, o~i are Pauli 
matrices, and e = u> corresponds to resonant excitation. 
The state vector is written as: 

We perform a rotating wave transformation by operating 
on with the unitary operator Q, where: 

Q = (cr + o*)/2 + exp(iwi + i(f>){<7 - a z )/2. (3) 

The Schrddingcr equation then takes the form (setting 

h = 1): |D = -iH(t)\i(t)} where the effective Hamilto- 
nian is given by: 

H = a{t)a+ + a*(t)a-, (4) 

with a(t) = —(go/2) [exp(— i2ut — i2<fi) + 1], and in the 
rotating frame the state vector is: 

\i( t )) = Q\m = ( gjjj ) . (5) 

Now, one may choose to make the rotating wave approxi- 
mation (RWA) , corresponding to dropping the fast oscil- 
lating term in a(t). This corresponds to ignoring effects 
(such as the Bloch-Siegert shift ) of the order of (go/oj), 
which can easily be observable in experiment if go is large 
jl|-||. On the other hand, by choosing go to be small 
enough, one can make the RWA for any value of lo. We 
explore both regimes in this paper. As such, we find the 
general results without the RWA. 

From EqsJJ and ||, one gets two coupled differential 
equations: 

C (t) = -(g /2) [1 + exp(-i2ut - i2<f>)] C x (t) (6a) 
Ci(t) = -(go/ 2) [l + ex.p(+i2wt + i2<t>)]C (t). (6b) 
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We assume |Co(t)| 2 = 1 is the initial condition, and pro- 
ceed further to find an approximate analytical solution of 
Eq.||. Given the periodic nature of the effective Hamilto- 
nian, the general solution to Eq|j| can be written in the 
form: 



l&)> = £ ( T )exp(nH2^-z20)). 



(J) 



Inserting Eq.0 in Eq.^|, and equating coefficients with 
same frequencies, one gets for all n : 



a n = i2nuja n + igo(b n + b n -i)/2, 
b n = i2nujb n + ig a {a n + a n+1 )/2. 



(8a) 
(8b) 



Here, the coupling between ao and bo is the conventional 
one present when the RWA is made. The couplings to 
the nearest neighbors, a±i and 6±i, are detuned by an 
amount 2u>, and so on. To the lowest order in (go /to), we 
can ignore terms with \n\ > I, thus yielding a truncated 
set of six equations: 



«o = igo(bo + b-i)/2, 

bo = igo{a Q + a x )/2, 

di = i2uja 1 + igo{b\ + b )/2, 

b± = i2ujbi + igoai/2, 
d_i = — z2wa_i + i<?ofe_i/2, 
&_i = — i2cij&_i + igo(a_i + Oq)/2. 



(9a) 
(9b) 
(9c) 
(9d) 
(9e) 
(9f) 



We consider go to have a time-dependence of the form 
go(t) = goM [1 — exp(— t/r sw )], where the switching time 
constant t sw is large compared to other characteristic 
time scales such as 1/lo and 1/ ' goM- Under this condition, 
one can solve these equations by employing the method 
of adiabatic elimination, which is valid to first order in 
r\ = (<7o/4w). Note that r\ is also a function of time, and 
can be expressed as r](t) = i]o [1 — exp(— t/r sw )], where 
Va = (9oM/4u). To solve the set of equations above, we 
consider first Eqs.^ejand 61 In order to simplify these two 
equations further, one needs to diagonalize the interac- 
tion between a_i and b-\. Define /i_ = (a_i — 6_i) and 
/x + = (a_i + b-i), which now can be used to re-express 
these two equations in a symmetric form as: 

A- = -i{2w + go/2)n- - ig a /2, (10a) 
M+ = -i(2u - g /2)n + + ig Q a /2. (10b) 

Adiabatic following then yields (again, to lowest order 
in rj ): fi- w — rjao and /i + ks rjao, which in turn yields 
a_i ~ and b—i ~ rjao- In the same manner, we can 
solve equations [k] and [k], yielding: a\ ~ — i]bo and 6i w 
0. 

Note that the amplitudes of a_i and b\ are vanish- 
ing (each proportional to rf ) to lowest order in rj, and 
thereby justifying our truncation of the infinite set of re- 
lations in Eq.^. It is easy to show now: 



ao = igobo/2 + iA(t)a /2, 
b = igoao/2 - iA(t)b /2, 



(11a) 
(lib) 



where A(t) = g^(t)/Au! is essentially the Bloch-Siegert 
shift. Eq.|ll| can be thought of as a two-level system ex- 
cited by a field detuned by A. For simplicity, we assume 
that this detuning is dynamically compensated for by ad- 
justing the driving frequency u>. This assumption does 
not affect the essence of the results to follow, since the 
resulting correction to r\ is negligible. With the initial 
condition of all the population in |0) at t = 0, the only 
non-vanishing (to lowest order in 77 ) terms in the solu- 
tion of Eq.[| are: 

a (t) w cos(g' (t)t/2), b {t) « i sm(g' (t)t/2), 
ai(t) w — i?7 sin(^ (*)*/ 2) , and w r}cos(g' (t)t/2), 

where 

g' (t) = l/tJ*go(t)dt =30 [1- (i/ r ™) _1 exp(-t/r sw )] . 

We have verified this solution via numerical integration 
of Eq.^| as shown later. Inserting this solution in EqJ|, 
and reversing the rotating wave transformation, we get 
the following expressions for the components of Eq.||: 

Co{t) = cos(g' (t)t/2) - 2nZsin{g' (t)t/2), (12a) 
d(t) - te-^ t+ ^[sm(g' Q (t)t/2) + 

+ 2r?S* cos(5o(*)*/ 2 )]) ( 12b ) 

where we have defined E = (i/2) exp(— i(2ut + 2(f))). 
To lowest order in 77, this solution is normalized at all 
times. Note that if one wants to carry this excitation on 
an ensemble of atoms using 7r/2 pulse and measure the 
population of the state |1) after the excitation terminates 
(at t — t when g'(r)r/2 = 7r/2 ), the result would be a 
output signal given by, 
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d( ffo (r), </>)| 2 = - [1 + 2 f ?siii(2wT + 2<j>)} 



(13) 



which contains information of both the amplitude and 
the phase of the driving field. This is our main result . 

A physical realization of this result can be appreciated 
best by considering an experimental arrangement of the 
type illustrated in Fig^. Here, a single group of atoms 
(e.g., atoms held in a dipole force trap) are subjected to 
a resonant, oscillating magnetic field, which is turned on 
adiabatically with a switching time-constant t sw , starting 
at t=0. After an interaction time of r , the population of 
the excited state (|1)) is determined instantaneously (i.e., 
with a time-constant much faster than 1 /u> and 1 / goM ) 
by coupling this state to an optically excited state (|2)) 
with a laser, and monitoring the resulting fluorescence. 
Such a measurement would correspond to the expression 
in Eq.|l3|. Defining the phase of the field at t = r to be 
6 T = lut + (j), the expression of Eq.^j| can be rewritten 
as: IC^r)! 2 = § [1 + 2r) sin(2^ T )] . 

In Fig.|^(a) we have shown the evolution of the excited 
state population \C\ (r) | 2 as a function of interaction time 
t, using the analytical expression of Eq.|l2|(b). 
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FIG. 1. Schematic illustration of an experimental arrange- 
ment for measuring the phase dependence of the population 
of the excited state (a) The microwave field couples the 
ground state (|0)) to the excited state (|1)). A third level, 
|2), which can be coupled to |1) optically, is used to mea- 
sure the population of |1) via fluorescence detection, (b) The 
microwave field is turned on adiabatically with a switching 
time-constant r sm , and the fluorescence is monitored after a 
total interaction time of r. 






FIG. 2. Illustration of the Bloch-Siegert Oscillation 
(BSO): (a) The population of state |I), as a function of the 
interaction time r , showing the BSO superimposed on the 
conventional Rabi oscillation, (b) The BSO oscillation (am- 
plified scale) by itself, produced by subtracting the Rabi os- 
cillation from the plot in (a), (c) The time-dependence of 
the Rabi frequency. Inset: BSO as a function of the absolute 
phase of the field. 

Under the RWA, this curve would represent the con- 
ventional Rabi oscillation. However, we notice here 
some additional oscillations, which is magnified and 
shown separately in figure 2(b), produced by subtract- 
ing the conventional Rabi oscillation (sin 2 (g(i)/2)) from 
figure 2(a). That is, figure 2(b) corresponds to what 
we call the Bloch-Siegert Oscillation (BSO), given by 
77 sin(gg (r)r) sin(2</> r ). The dashed curve (c) shows the 



time-dependence of the Rabi frequency. These analyt- 
ical results agree very closely to the results which are 
obtained via direct numerical integration of Eq|| Note 
that the BSO is at twice the frequency of the driving 
field, and its amplitude is enveloped by a function that 
vanishes when all the atoms are in a single state. 

Consider next a situation where the interaction time,r, 
is fixed so that we are at the peak of the BSO envelope 
(which corresponds to a (7r/2) pulse for the Rabi oscilla- 
tion). We further assume that r is long enough so that 
<7o(t) « goM- The experiment is now repeated many 
times, with a different value of (f> each time. The cor- 
responding population of |1) is given by 77 O sin(20 T ) = 
770 sin(2(cjT -I- </>)), and is plotted as a function of <j) in 
the inset of figure |[ This dependence of the population 
of |1) on the initial phase cf> (and, therefore, on the final 
phase 4> T ) makes it possible to measure these quantities. 

As indicated above, these effects can be observed in an 
experiment where, for example, a stationary collection 
of atoms are excited repeatedly by a microwave field. 
However, a more robust process for observing this effect 
can be realized using an atomic beam. For illustration, 
consider first a situation where the atoms are emitted 
in regular intervals At, and propagate in the z direction 
with a fixed velocity v. We can describe such a source by 
the line density of the number of atoms at a position z 
and at an instant t: M(z, t) = m 8[z — v(t — I At)] . 

The microwave field is assumed to be spatially vary- 
ing, corresponding to a Rabi frequency that vanishes 
for z < zq. For z > zq, it is given by: go(z) — 
9om [1 - exp(-(z - z )/z sw )] where z sw = vt sw repre- 
sents the distance over which the field is switched on. 
Under these assumptions, and in the limit where At — > 
(corresponding to a continuous, monovelocity atomic 
beam), it is easy to show that the normalized popu- 
lation S of |1) measured at a position z = z Q + vt 
(where r is the interaction duration corresponding to 
a 7r/2 pulse) as a function of time is given simply by 
[1/2 + 770 sin(2o;t + 2(f))], assuming that the microwave 
field at this position is B = Bocos(ujt + <j>). Note that 
(u)t + <f>) represents the absolute phase of the field as seen 
by the atoms when they arrive at this position. Thus, 
measurement of S directly reveals the absolute phase 
of the field, modulo 27r. Alternatively, one can mix S 
with another signal F corresponding to the second har- 
monic of a phase-shifted version of the microwave field 
(F = Fq cos[2(wi + <f> — 7r/2 — 8)]), where 9 is a controlled, 
variable phase shift), and observe the dc component of 
this signal, which will be proportional to cos(0). 

Consider next the more easily realizable situation 
where the atomic beam is produced by an effusive oven. 
Such a beam is typically characterized by a normalized 
velocity distribution f(v) = 2v 3 u~ 4 exp(— v 2 /u 2 ), where 
the u = ^j2kT /m is the most probable valocity, K is the 
Boltzman constant, T is the temperature and m is the 
mass of the atom [pi. The atoms that contribute to S 
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come from all different velocity groups. As such, each 
group experiences a different initial phase (i.e., phase at 
z — zq), and one might think that this would cause the 
signal to wash out. However, notice that S corresponds 
not to the initial phase, but rather to the phase at the 
observation point. Since the atoms contributing to S all 
have, by definition, arrived at this point at the same time, 
the signal will have the same time dependence, indepen- 
dent of the velocity group. Explicitly, the expression for 
S now becomes: 



S(t) = / dvf(v)[\sm(g' (T v )r v /2)\ 2 + 

Jv=0 

rj sin(g' (T v )T v ) sin(2w(t + (/>))]. 



(14) 



where t v = tu/v is the effective interaction time for the 
atoms with velocity v, and is a constant, corresponding 
to the interaction time for the most pro bable velocity u. 
The terms inside the integral in Eq.[14j are independent 
of t and 4> ' an d the effect of this averaging is simply to 
reduce the amplitude of the oscillatory signal observed. 

When a quantum bit (qubit) is represented by two non- 
degenerate states of a massive particle, it is necessary to 
apply a field at a frequency matching the energy differ- 
ence between these states, in order to produce an arbi- 
trary rotation of the qubit. In order to minimize the 
decoherence rate of such a qubit, one often chooses to 
use low energy spin transitions. In general, one is inter- 
ested in performing these transitions as fast as possible 
. As such, it is desirable to use a strong Rabi frequency. 
The ratio of the Rabi frequency to the qubit transition 
frequency is therefore not necessarily very small. One 
example of such a situation is already seen to occur in 
qubits based on trapped ions, for example [p|-|lO|]. Un- 
der this condition, one can not ignore the effect of the 
counter-rotating term. However, as we have shown here, 
the degree of excitation (e.g., the amplitude of the ex- 
cited state) depends not only on the product of the Rabi- 
frequency and the duration of the excitation, but also on 
the phase of the field at the time the interaction stops. 
Thus, one has to keep track of the phase of the excita- 
tion field at the location of the qubit [ pl]JT^ ] . In principle, 
the phase-tracking approach embodied in this paper it- 
self can be used for this purpose. Alternatively, one has 
to limit the strength of the Rabi frequency to a level dic- 
tated by the precision required of the particular qubit 
operation involved. We note that this effect is present 
for both direct radio-frequency excitation, as well as for 
indirect Raman excitation, which are functionally equiv- 
alent [p^|-|l5|. Finally, we point out that by making use 
of distant entanglement, this mechanism may enable tele- 
portation of the phase of a field that is encoded in the 
atomic state amplitude, for potential applications to re- 
mote frequency locking Jlp| |l9[| . 

In conclusion, we have shown that when a two-level 
atomic system is driven by a strong periodic field, the 



Rabi oscillation is accompanied by another oscillation at 
twice the transition frequency, and this oscillation carries 
the information about the absolute phase of the driving 
field. One can detect this phase by simply measuring only 
the population of the excited state. We have also shown 
how this effect may be observed using an atomic beam 
even if it has a substantial velocity spread. Finally, we 
have shown how this effect has to be taken into account 
in qubit operations. 
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